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The  influence  of  electronic  structure  and  doping  level  on  the  thermoelectric  transport  properties  of  BiSb3_ 
Te6  compound  was  investigated  using  first-principle  calculations  and  Boltzmann  transport  theory.  The 
calculation  results  indicate  that  BiSb3Te6  is  a  narrow  gap  semiconductor  with  an  indirect  band  gap  of 
0.1 13  eV.  The  calculated  band  structures  of  the  compound  present  nonparabolic  curves  with  multi-valley 
band  property.  The  density  of  states  is  enhanced  below  a  few  tenths  of  Fermi  energy  and  the  bands  are 
flatted  with  heavy  effective  mass.  The  calculated  thermoelectric  coefficients  of  BiSb3Te6  are  anisotropic 
and  mainly  dependent  on  the  doping  concentrations.  The  dopant  dependences  of  ZT  values  are  given  with 
the  optimal  carrier  concentrations  ~1019  cm-3.  This  theoretical  investigation  gives  a  valuable  insight  into 
the  relationship  between  the  electronic  structure  and  thermoelectric  transport  properties  of  BiSb3Te6 
material,  and  provides  a  quantificational  doping  level  to  improve  the  thermoelectric  performance. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  material  (TM)  provides  a  useful  solid-state  con¬ 
version  technology  for  refrigeration  and  power  generation  [1-4]. 
The  conversion  efficiency  of  TM  devices  is  limited  by  a  dimension¬ 
less  figure  of  merit  ZT  =  S2aTlk ,  where  S  is  the  Seebeck  coefficient,  T 
is  the  absolute  temperature,  a  is  the  electrical  conductivity,  and  k 
is  the  thermal  conductivity  which  generally  includes  the  contribu¬ 
tions  of  phonons  /q  and  electrons  ke.  High  performance  TM  requires 
a  large  power  factor  ( S2o )  and  a  low  thermal  conductivity.  Finding 
a  superior  ZT  material  is  full  of  challenges  due  to  the  conflicting 
thermoelectric  parameters  (S,  cr,  and  ke).  These  thermoelectric 
coefficients  are  highly  dependent  on  interrelated  material  carrier 
concentration  and  need  to  be  optimized  to  maximize  ZT  [5,6]. 

At  room  temperature  (up  to  500  I<),  BiSb3Te6  is  one  of  the 
best  performances  TMs  with  peak  ZT  values  of  0.96-1.4  [7,8]. 
To  further  improve  the  thermoelectric  property  of  BiSb3Te6  com¬ 
pound,  one  way  is  to  reduce  the  thermal  conductivity  of  the 
compound.  It  was  reported  that  a  minimal  /q  can  reach 
0.22Wm“1lr1  [9].  However,  the  obtained  ZT  was  less  than 
1.5  due  to  its  lower  power  factor  (PF)  of  inappropriate  carrier 
concentration.  Another  method  is  developed  to  optimize  the 
power  factor  by  adjusting  the  carrier  concentration.  A  maximum 
PF  of  7.21  pWcm-1  k-2  has  been  obtained  in  BiSb3Te6  material 
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[10].  In  addition,  a  small  quantity  of  Mn  [11],  Cu  [12]  and  Ag 
[13]  elements  were  introduced  into  BiSb3Te6  compound  to  mod¬ 
ulate  the  carrier  concentration.  Although  these  reports  suggested 
that  the  proper  doping  concentration  should  be  in  the  scope  of 
1019-1021  cm-3,  they  have  not  provided  the  carrier  concentration 
of  the  higher  power  factor  required.  The  quantitative  doping 
range  of  higher  ZT  is  still  lacking,  to  the  best  of  authors’ 
knowledgement. 

In  view  of  this,  the  purpose  of  this  study  is  to  provide  predic¬ 
tions  of  the  dependence  of  the  thermoelectric  coefficients  on  car¬ 
rier  concentration  with  the  Boltzmann  transport  theory  using  the 
full  electronic  structure,  giving  reasonable  doping  levels  to  opti¬ 
mize  the  thermoelectric  properties  of  BiSb3Te6  material. 

2.  Calculation  methods  and  theory 

2.1.  Calculated  methods 

The  crystal  structure  of  bulk  BiSb3Te6  was  constructed  substi¬ 
tuting  one  Bi  atom  for  Sb  atom  in  a  1  x  1  x  2  supercell  of  the  nom¬ 
inal  formula  Sb4Te6,  as  shown  in  Fig.  1.  To  confirm  the  stability  of 
the  crystal,  a  negative  formation  heat  of  -1.42  eV  was  obtained  by 
calculating  F(Bisb3Te6)  +  E(sb)~  2F(sb2Te3)  -  E(Bi).  The  optimized  lat¬ 
tice  constants  are  a  =  4.33  A  and  c  =  30.40  A,  which  are  very  close 
to  the  experimental  data  of  a  =  4.28  A  and  c  =  30.49  A  [14].  These 
results  suggest  that  the  BiSb3Te6  crystal  is  a  stable  compound 
and  can  be  synthesized  experimentally. 
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Fig.  1.  The  layered  (Bi,  Sb)4Te6  crystal  structure  (a)  was  constructed  from  the  Sb2Te3  unit  cell,  (b)  by  replacing  Sb  with  Bi  atom.  The  building  Brillouin  zone  (c)  and  the  high- 
symmetry  points  were  based  on  the  1x1x2  supercell. 


Based  on  the  density  functional  theory,  our  calculations  were 
performed  using  full-potential  linearized  augmented  plane-wave 
(LAPW)  [15]  method  as  implemented  in  WIEN2K  code  [16].  The 
generalized  gradient  approximation  with  the  Perdew-Burke-Ern- 
zerho  (PBE)  [17]  and  Engel-Vosko  (EV)  [18]  was  applied  to  evalu¬ 
ate  exchange-correlation  energy.  The  PBE-GGA  was  used  to 
calculate  the  structural  properties.  The  band-structure  properties 
were  computed  using  the  EV-GGA,  which  was  a  very  accurate 
function  in  determining  the  electronic  structure  of  thermoelectric 
materials  [19-21].  LAPW  spheres  of  radii  2.5  Bohr  were  employed 
for  both  Bi,  Sb  and  Te  atoms,  with  well-converged  basis  sets  deter¬ 
mined  by  min(ftMT)max(kn)  =  7.0,  where  RMT  is  the  minimum  LAPW 
radius  and  kn  is  the  plane-wave  cutoff.  Because  Bi,  Sb  and  Te  are 
heavy  elements,  the  relativistic  pi/2  local  orbital  [15]  is  introduced 
in  the  spin-orbit  (SO)  calculations. 

An  8  x  8  x  6  Monkhorst-pack  k-point  sampling  was  applied  in 
the  irreducible  Brillouin  zone  (BZ).  The  crystal  structure  was 
relaxed  by  the  BFGS  method  [22]  for  both  lattice  parameters  and 
atomic  internal  coordinates  until  the  stress  along  each  lattice 
direction  was  less  than  0.02  kbar,  and  the  total  energy  less  than 
0.0001  eV.  The  electronic  structure  was  calculated  in  the  optimized 
crystal  structure  with  32  x  32  x  24  denser  grids.  Thermoelectric 
transport  coefficients  were  derived  from  the  electronic  structure 
with  the  Boltzmann  transport  theory. 


2.2.  Boltzmann  transport  theory 

The  thermoelectric  parameters  were  calculated  by  resolving  the 
Boltzmann  transport  equation  with  the  constant  relaxation  time 
approximation  (CRTA)  [20].  The  transport  coefficients  of  electrical 


conductivity  cr,  Seebeck  coefficient  S,  and  electronic  thermal  con¬ 
ductivity  ke  are  given  by  the  following  equations  [19]: 


c(T,n)=Xm  (1) 

S(T,/|)=1(2X<°>)_1X«  (2) 

ke(T,n)=±tX^  (3) 

X(n)  =  1  J  Me)(e  -  (4) 


where  the  symbol  e  is  the  electrical  charge,  T  is  the  Kelvin  temper¬ 
ature,  Q  is  the  reciprocal  space  volume,  s  is  the  energy  to  the  corre¬ 
sponding  k  point,  p  is  the  chemical  potential,  and  /  is  the  Fermi 
distribution  function.  The  conductivity  tensors  [23]  is  written 
as: 

M£)  =  ^-^ya(k)Vf(k)  <?(£d£S|:)  (5) 

i,k 

where  t  is  the  relaxation  time,  va(k)  is  the  electrical  group  velocity, 
and  N  is  normalization  depending  on  the  number  of  k  points  sam¬ 
pled  in  the  BZ.  For  resolving  these  equations,  we  adopt  t  as  a  con¬ 
stant  without  variation  on  the  scale  of  a  few  kT  in  the  inelastic 
scattering  regime.  This  is  usually  applied  to  degenerately  doped 
semiconductors  and  metals,  obtaining  many  reasonable  results 
between  the  calculations  and  experiments  [24-28]. 


S.  Feng  et  al.  /  Computational  Materials  Science  95  (2014)  563-567 


565 


Fig.  2.  Electronic  DOS  and  band  structure  of  the  BiSb3Te6  crystal  calculated  using 
the  EV-GGA  along  the  high-symmetry  points  in  the  first  Brillouin  zone. 

3.  Results  and  discussions 

3.1.  The  electronic  structure 

The  calculated  density  of  states  and  band  structures  of  the 
BiSb3Te6  crystal  along  the  high  symmetry  lines  are  plotted  in 
Fig.  2.  It  indicates  that  the  crystal  is  a  narrow-gap  semiconductor 
with  an  indirect  band  gap  Eg  of  0.113  eV.  This  result  agrees  well 
with  the  experimental  data  of  0.112  eV  [29],  but  is  lower  than 
the  value  of  0.19  eV  30,31].  The  Fermi  level  approaches  the 
valence  bands,  with  the  implication  that  the  crystal  is  a  p- type 
semiconductor  with  the  hole-transport  property.  From  the  band 
structures  of  Fig.  2,  the  conduction-band  maximum  is  along  G-Z 
direction.  Along  this  direction,  a  relative  flat  band  is  observed, 
leading  to  the  formation  of  heavy  effective  mass  and  lower  carrier 
mobility.  Unlike  the  conduction-band  structures,  the  valence-band 
maximum  is  along  G-M  direction  and  the  Fermi  energy  is  close  to 
the  peak  point  of  the  valence  band.  With  dropping  the  energy, 
about  0.2  eV,  into  the  band  edge  of  G  point,  the  valence  bands 
become  heavier.  The  heavy  bands  are  attributed  to  the  strong 
effect  of  spin-orbit  interactions  in  the  BiSb3Te6  crystal,  which  is 
very  similar  to  the  Bi2Te3  compound  32,33]. 

In  addition,  the  density  of  states  (DOS)  rapidly  increases  with 
lowering  a  few  tenths  of  the  Fermi  energy  into  the  valence  bands. 
This  rapid  change  of  DOS  is  beneficial  to  improve  the  thermoelec¬ 
tric  power  of  the  TMs  [34,35].  Meanwhile,  the  BiSb3Te6  crystal  pre¬ 
sents  multivalley  band  structures  and  these  bands  are  highly 
nonparabolic,  which  are  contributable  to  the  enhancement  of  ther¬ 
moelectric  performance  at  elevated  temperature  and  doping  level 
[7,36]. 


Fig.  4.  Calculated  power  factor  ( PF )  as  a  function  of  chemical  potentials  for  both  n- 
type  and  p-type  doping  at  room  temperature.  The  reasonable  doping  levels  are  in  a 
region  between  the  half  values  of  the  maximum  PF  at  p1  and  n i  points. 

3.2.  The  transport  coefficients 

Fig.  3a  shows  the  calculated  room-temperature  Seebeck  coeffi¬ 
cients  in-layer  (along  the  cleavage  plane)  Sxx ,  and  cross-layer  (zz) 
Szz  with  both  n-type  and  p- type  doping.  The  Seebeck  coefficients 
of  p- type  doped  crystal  increase  rapidly  than  that  of  n-type  doped 
crystal,  reaching  the  maximum  values  of  253  pV  K_1  and 
-185  pVK-1.  They  are  very  close  to  the  experimental  results  of 
256  pVIC1  [6]  and  255  pVK”1  [13],  but  larger  than  that  of 
-149  pV  K_1  [6].  According  to  Mott  theory  [35],  the  improvement 
of  Seebeck  coefficients  closely  links  with  the  variation  of  the  DOS 
near  Fermi  energy.  Therefore,  higher  Seebeck  coefficients  in  this 
calculation  originate  from  the  rapidly  changing  DOS,  as  plotted  in 
Fig.  2.  This  enhancement  of  Seebeck  coefficients  has  also  been 
experimentally  proved  in  the  PbTe.Tl  compound  [37]. 

As  indicated  in  Fig.  3a,  we  also  observe  the  BiSb3Te6  crystal  dis¬ 
playing  isotropic  Seebeck  effect.  The  extreme  points  of  Seebeck 
coefficient  are  within  the  scope  of  the  energy  gap;  it  implies  that 
the  Seebeck  coefficient  can  be  improved  by  adjusting  the  Fermi 
energy.  In  contrast  to  the  Seebeck  coefficient,  the  electrical  conduc¬ 
tivity  is  anisotropic  for  both  p- type  and  n-type  materials.  The  ratio 
of  (Txxl(jzz  reaches  ~4.5  for  n-type  doping,  and  ~1.5  for  p-type  dop¬ 
ing,  as  shown  in  Fig.  3b.  Such  a  strong  anisotropic  electrical  con¬ 
ductivity  is  very  similar  to  the  layered  Bi2Te3  crystal  [38]. 

Fig.  4  is  the  comprehensive  transport  parameter  PF  with  the 
extreme  points  (pi,  p2,  n^,  n2)  located  at  different  chemical  poten¬ 
tial  for  both  p-  and  n-type  materials.  The  chemical  potentials  with 


Kl(eV)  |T(eV) 

Fig.  3.  Calculated  room-temperature  transport  coefficients  of  Seebeck  coefficient  (a)  and  electrical  conductivity  and  (b)  are  plotted  as  a  function  of  the  chemical  potential  for 
both  n-type  and  p-type  doping.  Symbols  of  xx  and  zz  are  parallel  and  perpendicular  to  the  cleavage  plane,  respectively. 
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Fig.  5.  The  dopant  dependence  of  the  room-temperature  thermoelectric  coefficients:  S,  o  and  PF  (a),  ZT  and  electronic  ke  (b). 


respect  to  pi  and  nA  points  are  -0.101  eV  and  0.074  eV,  respec¬ 
tively.  These  two  peak  points  are  within  the  feasible  doping  levels 
because  they  approach  the  edge  of  the  conduction  band  (n-type) 
and  valence  band  (p-type).  However,  p2  and  n2  points  are  located 
at  heavy  doping  concentrations  (~1021  cm-3),  in  which  a  large  ke 
is  introduced  into  TMs  and  thus  lowers  the  ZT  values.  In  this  calcu¬ 
lation,  the  available  doping  levels  therefore  are  in  a  region  between 
the  half  values  of  the  maximum  PF  at  pi  and  ni  points.  They  can  be 
limited  to  a  range  between  -0.15  eV(p- type)  and  +0.12  eV(n-type), 
corresponding  to  the  carrier  concentration  from  +2.8  x  102Ocm-3 
to  -8.5  x  1019  cm-3. 


Additionally,  the  average  power  factors  using  the  data  calcu¬ 
lated  from  the  xx-  and  zz-planes  are  drawn  with  the  solid  line  in 
Fig.  4.  The  optimized  PF  at  p 1  point  reaches  2.25  x  1011  W  m-1  - 
K-2  s-1,  about  7.2  x  10-3  W  m-1  K-2  if  t  is  set  at  3.2  x  10-14  s,  in 
good  agreement  with  the  experiment  data  of  7.24  W  m-1 1<-2  by 
Zhao  et  al.  [10].  Introducing  the  experimental  thermal  conductivity 
of  1.6  Wm-1  K-1  [39,40],  at  300  K,  the  average  dimensionless  fig¬ 
ure  of  merit  ZT  could  reach  1.35,  closes  to  1.2  as  reported  by  Poudel 
et  al.  41  ].  At  this  optimal  doping  concentration,  the  ZT  value  could 
be  further  improved  in  the  even  lower  thermal-conductivity  nano¬ 
scale  compound. 


Fig.  6.  Calculated  in-plane  electrical  conductivity  (a),  Seebeck  coefficient  (b)  and  power  factor  (c)  as  a  function  of  carrier  concentration  in  p-type  doped  BiSb3Te6  materials  at 
300-500  K. 
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3.3.  The  thermoelectric  properties  at  different  doping  levels 

Fig.  5a  displays  the  carrier  concentration  dependence  of  thermo¬ 
electric  parameters  for  the  p- type  doped  compound  at  room  tem¬ 
perature.  It  implies  that  with  increasing  the  doping  concentration, 
the  electrical  conductivity  continuously  increases  and  the  Seebeck 
coefficient  first  rises  and  then  reduces.  As  the  doping  concentration 
is  over  9.8  x  1018  cm-3,  the  Seebeck  coefficient  and  electrical  con¬ 
ductivity  change  in  opposite  direction  with  increasing  the  doping 
levels,  resulting  that  the  optimum  power  factor  does  not  match 
with  the  optimum  Seebeck  coefficients.  Here,  the  maximum  PF 
can  reach  2.78  x  1011  W  m-1 1<-2  s-1  at  n  =  1.06  x  1020  cm-3. 

High  carrier  concentration  BiSb3Te6  semiconductors  have  large 
power  factor.  However,  high  carrier  concentration  also  results  in 
large  electronic  thermal  conductivity  ke.  As  shown  in  Fig.  5b,  the 
electronic  thermal  conductivity  ke  attains  0.4  W  m_1  K1  at  the 
doping  level  of  5.5  x  1019cm-3,  approaching  the  same  order  of 
magnitude  of  the  lattice  thermal  conductivity  0.7  W  m-1  K1 
[30].  By  combining  the  available  experimental  thermal  conductiv¬ 
ity  42]  with  our  calculated  PF,  the  largest  ZT  value  is  up  to  1.64  at 
the  carrier  concentration  of  3.5  x  1019  cm-3.  This  predicted  figure 
of  merit  (ZT)  is  in  good  agreement  with  the  experimental  result 
of  1.56  reported  by  Jiang  et  al.  [8].  If  the  ZT  value  need  exceed 
1.5,  the  optimized  carrier  concentration  should  be  within  the  dop¬ 
ing  range  2.32  x  1019-4.55  x  1019  cm-3.  Additionally,  assuming 
that  the  lattice  thermal  conductivity  /q  may  suppress  small  enough 
(~0.8  W  m-1  K1)  at  room  temperature,  the  ZT  value  of  the  p- type 
doped  BiSb3Te6  compound  can  break  through  3,  larger  than  the 
currently  reported  ZT  =  2.6  in  SnSe  crystals  [43]. 

In  general,  the  doped  BiSb3Te6  materials  are  primarily  used  for 
thermoelectric  generator  at  300-500  K  [44].  Fig.  6  illustrates  the 
effect  of  doping  levels  on  the  thermoelectric  transport  properties 
of  the  doped  BiSb3Te6  materials  along  the  cleavage  plane  at  elevated 
temperatures.  It  can  be  observed  that  the  temperatures  have  an 
influence  on  the  electrical  conductivity  and  Seebeck  coefficient  at 
the  doping  level  below  1.5  x  1019cm-3.  After  that,  as  shown  in 
Fig.  6a,  the  electrical  conductivity  is  independent  of  the  temperature 
and  grows  up  linearly  with  increasing  the  carrier  concentration. 

In  heavy  doping  region  (>102°  cm-3),  as  shown  in  Fig.  6b,  the 
increase  in  temperature  results  in  a  high  Seebeck  coefficient.  The 
accelerated  rate  of  Seebeck  coefficient  with  temperature  can  reach 
0.2  qV  I<“2.  Since  the  PF  ( S2a )  is  more  sensitive  to  the  Seebeck  coef¬ 
ficient  than  electrical  conductivity,  as  displayed  in  the  dot  lines  of 
Fig.  6c,  the  extreme  point  of  the  PF  lifts  to  the  high  doping  level 
with  increasing  the  temperature.  The  optimized  PF  values  are 
therefore  improved  with  the  elevating  temperature  and  carrier 
concentration.  This  variation  of  the  optimal  PF  yields  functionally 
graded  thermoelectric  materials  that  are  available  for  improving 
the  coefficient  of  performance  of  the  thermoelectric  power 
generator. 


4.  Conclusions 

First-principle  calculations  were  used  to  investigate  the  elec¬ 
tronic  structure  and  thermoelectric  properties  for  the  BiSb3Te6 
alloy.  Similar  to  Bi2Te3,  BiSb3Te6  is  a  narrow-gap  semiconductor 
with  an  indirect  band  gap  of  0.1 13  eV,  about  0.05  smaller  than  that 
of  Bi2Te3.  The  band  structures  present  multivalley  characteristics. 
The  Fermi  energy  is  close  to  the  valence-band  maximum  and  with 
lowering  energy  ~0.2  eV;  the  density  of  states  for  the  crystal 
improves  greatly  and  the  valence  bands  become  flat.  Calculated 
thermoelectric  coefficients  of  the  BiSb3Te6  material  perform  aniso¬ 
tropic  properties.  Along  the  cleavage  plane,  maximum  Seebeck  coef¬ 
ficients  of  253  pV  K1  and  -185  pV  K-1  are  obtained  for  both  p-type 
and  n-type  materials,  respectively.  The  electrical  conductivity  of  the 


material  in-layer  is  ~4.5  times  larger  than  that  of  cross-layer  for  the 
n-type  doping,  and  ~1.5  times  for  p-type  doping.  To  acquire  higher 
performance  of  doping  BiSb3Te6  materials,  the  optimized  carrier 
concentration  should  be  within  the  range  2.32  x  1019  to 
4.55  x  1019  cm-3. 
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